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ABSTRACT 

A linear modulation of the primordial perturbations is proposed as an ex- 
planation for the observed asymmetry between the northern and southern hemi- 
spheres of the Wilkinson Microwave Anisotropy Probe ( WMAP) data. A cut 
sky, reduced resolution third year "Internal Linear Combination" (ILC) map 
was used to estimate the modulation parameters. A foreground template and 
a modulated plus unmodulated monopole and dipole were projected out of the 
likelihood. The effective chi squared was reduced by nine for three extra pa- 
rameters. The mean galactic colatitude and longitude, of the modulation, with 
68%, 95% and 99.7% confidence intervals were 56±i?±||±ff and 63±1^±^. The 
mean percentage change of the variance, across the pole's of the modulation, was 
62tjgi35t47- Implications of these results and possible generating mechanisms 
are discussed. 

Subject headings: cosmic microwave background, early universe, large-scale struc- 
ture of universe 



1. Introduction 

A fundamental assumption of cosmology is that the Universe is isotropic. This was 
confirmed, for the mean temperature of the cosmic microwave background (CMB), by the 
FIRAS experiment on the COBE satellite (Wright et al. 1992; Bennett et al. 1996). How- 
ever, the higher precision, measurements from the Wilkinson Microwave Anisotropy Probe 
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( WMAP) satellite (Bennett et al. 2003; Hinshaw et al. 2006; Jarosik et al. 2006; Page et al. 
2006; Spergel et al. 2006), have an anomalously asymmetric distribution, in the temperature 
fluctuation statistics, between the northern and southern hemispheres of the sky (Eriksen 
et al. 2004; Hansen et al. 2004a; Vielva et al. 2004; Park 2004; Copi et al. 2004; Hansen 
et al. 2004b; Larson & Wandelt 2004; Cruz et al. 2005; Land & Magueijo 2005a; Hansen 
2004; Bernui et al. 2005, 2006). On scales greater than about 5°, the variance of the CMB 
temperature fluctuations is anomalously higher in the southern hemisphere, in both galactic 
and ecliptic coordinates, compared to the northern hemisphere (Eriksen et al. 2004; Hansen 
et al. 2004a). This asymmetry also appears in higher order statistics (Vielva et al. 2004; 
Park 2004; Copi et al. 2004; Hansen et al. 2004b; Larson & Wandelt 2004; Cruz et al. 2005; 
Land & Magueijo 2005a; Hansen 2004; Bernui et al. 2005, 2006). 

In a spherical harmonic representation, scales ranging from £ = 2 to £ = 40 were found 
to be asymmetric. When optimized over direction, only 0.3% of isotropic simulations were 
found to produce higher levels of asymmetry (Hansen et al. 2004a). 

The result is not sensitive to the frequency band of the CMB (Hansen et al. 2004a) and 
a similar pattern (at lower significance) is seen in COBE (Hansen et al. 2004a). This argues 
against a foreground or systematics explanation. 

Although a simple single field inflation model would give isotropically distributed per- 
turbations, this is not necessarily the case in multi- field models (Linde & Mukhanov 2006). 
Thus, if it can be shown that the CMB fluctuations are not isotropic, it may be an indication 
that inflation was a multi-field process. 

The layout of the paper is as follows: In Sec. 2, a linearly modulated primordial power 
spectrum is proposed as the source of the observed isotropy breaking. Then, in Sec. 3, a 
method of evaluating the linear modulation parameters is outlined. The constraints are given 
in Sec. 4 and their implications and relation to other results are discussed in Sec. 5. 



2. Linear Modulation 

An isotropy-breaking mechanism may be parameterized as (Prunet et al. 2005; Gordon 
et al. 2005; Spergel et al. 2006) 

<5T(n)=5T iso (n)(l-/(n)) (1) 

where 5T is the observed CMB temperature perturbations, 5T iso are the underlying isotrop- 
ically distributed temperature perturbations and n is the direction of observation. An 
isotropic distribution of perturbations is recovered when / = 0. 
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Spergel et al. (2006) parameterized / as 

i i 

=EE f*nYU*) (2) 



=1 m= 



where Y^ m are spherical harmonics and j = 1 and j = 2 were tried. As / will be a real 
function, the condition f^ m = {—l) m fl_ m is required, where the asterisk indicates the com- 
plex conjugate. The effective chi squared improvement, A-xls = — 2Alog£ is only -3 for the 
j = 1 case and only -8 for the j = 2 case, where C is the likelihood (Spergel et al. 2006). 

The dipolar modulation (j = 1) case has the potential to explain the lack of isotropy be- 
tween two hemispheres, as it will reduce the variance of the perturbations in one hemisphere 
and increase it in the other. 

In this article, an underlying spatial model for a dipolar modulation is investigated. On 
large scales (£ < 30) the main contribution (Hu & Sugiyama 1995) to the CMB perturbations 
is the Sachs Wolfe effect (Sachs & Wolfe 1967) 

8T(h) 

where $ is the curvature perturbation, in the Newtonian gauge (Bardeen 1980), evaluated 
at last scattering. 1 The left hand side is the measured temperature fluctuation divided by 
the average measured temperature. 

It follows that on large scales and at last scattering 

$ = (l-/)$ iso . (4) 

A dipolar modulation of the last scattering surface could result from a spatially linear mod- 
ulation 

/(n)=w-n (5) 

where w = (w x ,w y ,w z ) is the gradient of the modulating function and the dot indicates a 
dot product. 

The spherical harmonic representation is related to a Cartesian representation by 

[/ 10 , Re(/n), Im(/ n )] = ^| ^2w z , -V2w x , V2w y ~\ (6) 

As they are linearly related, a uniform prior on one translates into a uniform prior on the 
other. 



1 This is the surface where the Universe becomes effectively transparent to photons and occurs at a redshift 
of about 1000. 
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3. Likelihood Analysis 

The likelihood is assumed to follow a multivariate Gaussian distribution 

^oclCr^exp^TTC-V) (7) 

where T is a vector of the temperature measurements in unmasked areas of the CMB pix- 
elized map. Each element of the covariance matrix, C, is evaluated by 

C(n, m) = (1 - /(n))C iso (n, m)(l - /(m)) + AC marg (n, m) (8) 

where C(n, m) corresponds to the covariance between pixels at position n and m. The 
instrumental noise is negligible on large scales (Jarosik et al. 2006; Hinshaw et al. 2006) and 
so is not included. The underlying isotropic covariance matrix, between pixels in directions 
n and m, can be decomposed as 

C iS o(n,m) =J2Cew 2 e P e (h - ih)(2£ + . (9) 

£ 

where Ce is the angular power spectrum, Pe is the Legendre polynomial of order £ and we 
is the effective window function of the smoothed map evaluated at a low pixel resolution. 
Modulating the temperature perturbations leads to a transformed covariance matrix of the 
form 

C(n, m) = (1 - /(n))C iso (n, m)(l - /(&)) . (10) 

A marginalization term (Tegmark et al. 1998; Bond et al. 1998; Slosar et al. 2004; Slosar 
& Seljak 2004; Hinshaw et al. 2006) for foregrounds and a modulated and unmodulated 
monopole and dipole was also added 

Cmarg = T^^ + + 2^ _ ^)( < ^° + ~ f) + ^foregrounds • (H) 

The unmodulated monopole and dipole terms are needed to account for any residual effects 
of the background temperature and peculiar motion of the observer (Sachs & Wolfe 1967). 
When A is made sufficiently large, the likelihood becomes insensitive to any terms included 
in C marg (Tegmark et al. 1998; Bond et al. 1998; Slosar et al. 2004; Slosar & Seljak 2004; 
Hinshaw et al. 2006). 

Third year WMAP data were used 2 and the preprocessing followed was the same as in 
the WMAP analysis (Hinshaw et al. 2006; Spergel et al. 2006) for their large scale likelihood 



2 Obtained from http://lambda.gsfc.nasa.gov/product/map/current/ 
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evaluation. The A^ side = 512 "Internal Linear Combination" (ILC) map 3 was masked with 
the Kp2 mask and smoothed with a 7.3° (FWHM) Gaussian smoothing function. It was 
then degraded to N s ^ c = 8 using the HEALPix 4 software package (Gorski et al. 2005). The 
Kp2 mask, consisting of zeros and ones, was also degraded to N si d c = 8 and any pixels with 
values larger than 0.5 were set to one, else they were set to zero. The smoothed degraded 
ILC map was then remasked with the degraded Kp2 mask. The foreground template was 
taken to be the difference between the raw V-band map and the ILC map. The foreground 
marginalization term in Eq. (11) was set equal to the outer product of the foreground tem- 
plate with itself (Tegmark et al. 1998; Bond et al. 1998; Slosar et al. 2004; Slosar & Seljak 
2004; Hinshaw et al. 2006). 

The C e values were treated as free parameters for £ = 2 to 10. For £ = 11 to 32, the 
unmodulated maximum likelihood values were used. The smoothing and degrading of the 
data make the likelihood insensitive to £ > 32. 

The likelihood, for the modulated model, was numerically maximized using a quasi- 
Newton method. 

Marginalized distributions of the parameters were obtained using the Metropolis al- 
gorithm. After an initial burn in run, a proposal covariance matrix was constructed from 
20000 samples. This was used to generate three additional sets of 20000 samples, each with 
a different starting value chosen from a burned in chain. The Gelman and Rubin test were 
used to check convergence and then the 60000 samples were used to evaluate the marginal 
probability distributions of the parameters. All priors were taken to be uniform. The upper 
and lower limit of each confidence interval was chosen so as to exclude the same number of 
samples above and below the interval. 



4. Results 

The improvement, in the likelihood, compared to an unmodulated model was Axl s = —9 
with three extra parameters (u>i, u>2, ^3). A plot of the maximum likelihood modulation 
function /, is shown in Figure 1. The marginalized distributions of C2 to C\o were found 
not to be significantly different from those in an unmodulated model (Hinshaw et al. 2006). 
The weight vector samples (w x ,w y ,w z ) were converted into Galactic co-latitude (varying 
between 0°,180°), longitude (varying between —180° and 180°) and percentage change of 



3 number of pixels =12-/V s 2 idc . 
4 http://hcalpix.jpl. nasa.gov 
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Fig. 1. — Maximum likelihood linear modulation function projected on the last scattering surface. High 
values indicate suppression of the primordial perturbations and low values indicate enhancement. 
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temperature variance in the direction of symmetry breaking (A = 200 (vol + w l + w z))- The 
two dimensional confidence intervals for the co-latitude and longitude are shown in Figure 2. 
The marginalized one dimensional distributions are shown in Figure 3. Only 0.2% of the 
samples were more than 90° from the maximum likelihood point. The results are summarized 
in Table 1. 



5. Discussion and Conclusions 

In this article the modulation model investigated by Spergel et al. (2006) has been 
extended by including a marginalization over the unmodulated monopole and dipole. This 
additional feature is required if the apparent isotropy breaking had a primordial origin. 
Including this marginalization improved the Axl s value from -3 to -9. 

As seen from the confidence intervals in Table 1 and Figure 3, the marginalized posterior 
probability of A has its maximum more than three sigma away from the unmodulated case 
(A = 0). The modulated model is also preferred by the Akaike Information Criteria (AIC) 
(Akaike 1974; Magueijo & Sorkin 2006). It is not preferred by the Bayesian Information 
Criteria (BIC) (Schwarzy 1978; Magueijo & Sorkin 2006). However, the BIC is an approxi- 
mation of the Bayesian evidence and assumes a prior for the parameters which is equivalent 
to one observation (Raftery 1995). The Bayesian evidence will be inversely proportional 
to the volume of the prior probability distribution of the modulation parameters. It may 
be hard to produce a modulation larger than one without effecting the observed dipole. 
A reevaluation of the Bayesian evidence is needed to see how it depends on the assumed 
prior. This could be done using a nested sampling algorithm (Mukherjee et al. 2006) which, 
unlike the BIC, does not require a Gaussian approximation to be made for the posterior 



Table 1. Marginalized Statistics 



Parameter 


Mean a 




A 


62lj* 


-35+57 
-35-47 


57 


Colatitude 




-36+65 
-35-51 


51 


Longitude 


fi o+28+59+105 
uo -26-58-213 


68 



including confidence intervals (68%, 95%, 99. 7%). 
b Maximum Likelihood. 
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Fig. 2. — The contours enclosing 68% and 95% of the Monte Carlo samples. The maximum likelihood (o) 
Galactic co-latitude and longitude is (51°, 68°). One of the most asymmetric directions found by Hansen 
et al. (2004a) (A) is (80°, 57°). The north ecliptic (□) is at (60°, 96°). This is plotted over the third year 
iV S idc = 512 ILC map (Hinshaw et al. 2006). 




20 40 60 80 100 120 50 100 150 -100 100 

A Co - latitude Longitude 



Fig. 3. — Marginalized distributions of the percentage change in the variance across the poles (A) and the 
galactic coordinates of the direction of the modulation. The 68%, 95% and 99.7% confidence intervals are 
also shown. 
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distribution. 

Spergel et al. (2006) also evaluated whether there was an additional quadrupolar com- 
ponent to the modulation. This component could potentially be useful in explaining the 
alignment and planarity of the quadrupole {£ = 2) and octopole {£ = 3) seen in the WMAP 
temperature data (de Oliveira-Costa et al. 2004; Schwarz et al. 2004). The normal direction 
of the plane of alignment is (30°,— 100°). Also, when the coordinate system is rotated in 
the direction of the normal of the £ = 2,3 planarity there is anomalous power in the m = 3 
component of the t = 5 multipole (Land & Magueijo 2005b). Spergel et al. (2006) found 
that including a dipolar and quadrupolar component to the modulation improved Axlg by 
only 8 for a total of 8 extra parameters. Higher order terms in a spatial modulation could be 
implemented as terms quadratic in the spatial coordinates. Whether these additional terms 
will become significant when an unmodulated monopole and dipole are marginalized over 
will be part of a future investigation. 

The effect of marginalization over foregrounds was checked and found not to play a big 
role. Similar improvements are obtained when the foreground corrected V band is used in- 
stead of marginalization. Also, the results are not sensitive to the exact method of degrading 
and applying the mask. A Kp2 extended mask (Eriksen et al. 2006) did not make a signifi- 
cant difference. Including additional Ci, with £ > 10, as parameters to be estimated (rather 
than set to their unmodulated ML values), also does not significantly effect the results. 

It is interesting to compare the estimated modulation found in this article to that 
of Hansen et al. (2004a). The 10 most effective axes of symmetry breaking, for a range 
of scales, are plotted in their Figure 24. A similar area, to the two dimensional confidence 
intervals in Figure 2, is covered. Also, their Figure 19 compares the power spectra in different 
hemispheres. The range of values is consistent with the confidence intervals for A in Table 1 
and Figure 3. 

Prunet et al. (2005) tested for a dipolar modulation. However, the largest scale they 
looked at was £ = 20 to 100 binned. They did not get significant results in that range. As 
the observed modulation only occurs for £ < 40 (Hansen et al. 2004a), the £ = 20 to 100 
range would not be expected to show significant modulation when binned. 

Freeman et al. (2006) propose that the modulation of £ = 2 to £ = 7 may be sensitive to 
any residual unmodulated dipole component. This is not a concern for the approach taken 
in this article as an unmodulated dipole is projected out of the likelihood, see Eq. (11). 

Searches for lack of isotropy using a method based on a bipolar expansion of the two 
point correlation function do not detect the north south asymmetry in the £ = 2 to £ = 40 
range (Hajian & Souradeep 2006; Armendariz-Picon 2006). The linear modulation model 
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could be used to understand why the bipolar estimator is insensitive to this type of isotropy 
breaking. 

A small scale cut off in the modulation implies that a linear modulation of the primordial 
power spectrum would only apply to wave numbers larger than about 4 x 10~ 3 /i Mpc -1 . It 
would be interesting to evaluate whether this modulation would be detectable in future large 
scale galaxy surveys. However, at a redshift of one the change in the variance at opposite 
poles would only be about four percent, due to the smaller comoving distance. 

A number of attempts have been made to explain the asymmetry in terms of local 
nonlinear inhomogeneities (Moffat 2005; Tomita 2005a,b; Inoue & Silk 2006). It would be 
interesting to see if the polarization maps of the CMB (Page et al. 2006) could be used to 
distinguish local effects from a modulation of the primordial perturbations. 

Primordial magnetic fields (Durrer et al. 1998; Chen et al. 2004; Naselsky et al. 2004), 
global topology (de Oliveira-Costa et al. 2004; Kunz et al. 2006), and anisotropic expansion 
(Berera et al. 2004; Buniy et al. 2006; Gumrukcuoglu et al. 2006) can also lead to isotropy 
breaking. However, in these cases the modulating function is of higher order than dipolar 
and so these mechanisms are better suited for explaining the alignment between i = 2 and 
£ = 3 and the high (£,m) = (5,3) mode (Gordon et al. 2005). 

An additive template based on a Bianchi Vllh model has been shown to provide a good 
fit to the asymmetry (Jaffe et al. 2005). However, the model is only empirical as it would 
require a very open Universe which is in conflict with many other observations. It is harder 
for additive templates to explain the alignment between 1 = 1 and i = 3 as this requires 
a chance cancellation between an underlying Gaussian field and a deterministic template 
(Gordon et al. 2005; Land & Magueijo 2006). 

As seen in Figure 2, the maximum likelihood direction of modulation was found to be 
about 44° from the ecliptic north pole. Only about 9% of the time would two randomly chosen 
directions be as close, or closer, together. This may be an indication that the modulation 
is caused by some systematic effect or foreground. However, as can be seen from Figure 2, 
the confidence intervals, for the direction of modulation, cover just under half the northern 
hemisphere. Therefore, the actual direction, of modulation may be significantly further away 
from the ecliptic north pole. 

Standard single field inflation would produce isotropic perturbations. However, multi- 
field models, such as in the curvaton scenario (Lyth & Wands 2002; Mollerach 1990; Linde 
& Mukhanov 1997; Enqvist & Sloth 2002; Moroi & Takahashi 2001), can produce, what to 
a particular observer appear to be, non- isotropic perturbations (Linde & Mukhanov 2006). 
The curvaton mechanism produces a web like structure in which relatively stable domains 
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are separated by walls of large nonlinear fluctuations. If the mass of the curvaton field is 
sufficiently small, our observable Universe could be enclosed within a stable domain. If we 
happen to live near one of the walls, of a domain, then the amplitude of the perturbations 
will be larger on the side of the observed Universe closer to the wall (Linde & Mukhanov 
2006). However, if our observed Universe was far enough away from the web walls, the very 
large scale fluctuations would be linear and so isotropy would be unlikely to appear to be 
broken (Lyth 2006). As the non-isotropic nature only extends to about I = 40 (Hansen et al. 
2004a), it would be necessary for the inflaton perturbations to dominate over the curvaton 
ones for wave numbers larger than about 4 x 10~ 3 h Mpc -1 . The curvaton produces curvature 
perturbations proportional to U 1 / 2 (Lyth & Wands 2002), where V denotes the inflaton 
potential. While the inflaton produces curvature perturbations proportional to V 3 ^ 2 /V, 
where V denotes the slope of the potential. So if there is a sudden drop in V and V, it is 
possible for the non-isotropic curvaton perturbations to dominate for wave numbers smaller 
than 4 x 10~ 3 h Mpc -1 and inflaton perturbations to dominate for larger wave numbers. 

There are oscillations in the WMAP power spectrum, at around I = 40, which may be 
caused by a change of slope in the inflaton potential (Covi et al. 2006). Whether all these 
elements can be put together to make a working curvaton model, that fits the data as well 

linear modulation, is still being investigated. 

The results presented here provide a parameterization for the observed asymmetry be- 
tween different hemispheres of the WMAP data. Having a specific model for the primordial 
fluctuations will make it easier to develop new tests for this asymmetry and help determine 
if it is a genuine window into new physics at the largest observable scales. 

I thank Olivier Dore, Wayne Hu, Dragan Huterer, Nemanja Kaloper, David Lyth and 
Hiranya Peiris for useful discussions. I was supported by the Beecroft Institute for Particle 
Astrophysics and Cosmology and also by the Kavli Institute for Cosmological Physics under 
NSF PHY-01 14422. Some of the results in this paper have been derived using the HEALPix 
(Gorski et al. 2005) package. The WMAP data were obtained from the Legacy Archive for 
Microwave Background Data Analysis (LAMBDA). 
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